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In a heterogeneous wireless cellular network, each user may be covered by multiple access points such as 
macro/pico/relay/femto base stations (BS). An effective approach to maximize the sum utility (e.g., system through- 
put) in such a network is to jointly optimize users' linear procoders as well as their base station associations. In this 
paper we first show that this joint optimization problem is NP-hard and thus is difficult to solve to global optimality. 
C/5 ■ To find a locally optimal solution, we formulate the problem as a noncooperative game in which the users and 

the BSs both act as players. We introduce a set of new utility functions for the players and show that every Nash 
£NJ ' equilibrium (NE) of the resulting game is a stationary solution of the original sum utility maximization problem. 

, Moreover, we develop a best-response type algorithm that allows the players to distributedly reach a NE of the 
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game. Simulation results show that the proposed distributed algorithm can effectively relieve local BS congestion 
and simultaneously achieve high throughput and load balancing in a heterogeneous network. 

I. Introduction 

Consider a multicell heterogeneous network (HetNet) in which every cell is installed with not only a macro base 
station (BS) but also a set of micro/pico/femto base stations, all equipped with multiple antennas and sharing the 
same frequency band. In such a network, the users are often simultaneously covered by multiple BSs with different 
capabilities and load status. If the users in a HetNet are assigned to BSs simply according to their received signal 
strengths (as is done in conventional cellular networks), then a BS located close to a hot spot may experience severe 
congestion, causing poor quality of service in the network. Indeed, it has been shown that the signal strength based 
approach for user-BS association in a HetNet can be highly suboptimal for congestion management and fairness 
provisioning (2l, Q- For a HetNet, the overall system performance depends not only on the physical layer choices 
of precoder design and power control, but also on how the users and the BSs are associated. Consequently, BS 
assignment should be an integral part of physical layer optimization of the overall system performance. 
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Conventionally, physical layer resource management and network performance optimization involve such aspects 
as transceiver design, power control and spectrum management, while the user-BS assignment is assumed known 
and fixed. For example, under this assumption, game theoretic approaches have been used to design optimal transmit 
precoders for a MIMO interference channel (MIMO-IC). The authors of 0], 0, among others, propose to formulate 
the transmit co variance matrix optimization problem as a noncooperative game in which the transmitters/users 
compete with each other for transmission. Reference ||6] further takes into consideration the robustness issue of the 
problem. In these studies, each user selfishly maximizes its own transmission rate, while treating all other users' 
interference as noise. Simple distributed algorithms (based on iterative water-filling) are derived with convergence 
guarantees, but the resulting solutions of the game are inefficient in terms of system throughput. This is due to both 
the lack of coordination among the transmitters/users and the choice of utility functions which do not consider the 
interference effect on other users in the system. 

In addition to the competitive design, we can also design transmit precoders by maximizing a suitable system 
utility function. Unfortunately, most optimization problems in this category have been proven to be NP-hard in 
various settings iPTI- fTPl . As a result, many authors focus on developing efficient algorithms to compute high 
quality sub-optimal solutions for these problems, e.g., in MIMO-IC, 031, El, MISO-IC Ql-Ql and MIMO-IC 
with a single data stream per user ||T6l , |H71 . In particular, references ifTTI . ITT21 propose an iterative algorithm based 
on the first order Taylor approximation of the nonconcave part of the weighted sum rate (WSR) objective. It is 
shown that the WSR values generated by this algorithm increase monotonically, but the convergence of the users' 
transmit covariances (to a stationary solution) is left unproven. For the problem of maximizing a general utility 
function, reference |J8] proposes a cyclic ascent method for linear precoder design in a MISO network. The proposed 
algorithm can deal with any smooth utility functions and is known to converge to a stationary solution. Reference 
|9l develops an algorithm that optimizes the Max-Min utility in a MIMO network by adapting the transmit and 
receive precoders alternately. References 11181 . ||T9l propose a weighted Minimum Mean Square Error (WMMSE) 
algorithm in which the transmitters and receivers iteratively update their linear transmit and receive strategies to 
optimize the system utility function. The authors show that as long as the system utility function satisfies some 
mild regularity conditions, their algorithm is guaranteed to converge to a stationary point of the problem. 

All of the above cited works aim at optimizing the linear transceiver structures under the assumption that the 
transmitter-receiver association is known and fixed. The problem of joint cell site selection and power allocation 
in the traditional CDMA based network has been first considered in |[20l - |[22l and later in a game theoretical 
perspective in ||23l , 11241 . The objective of the network optimization is to minimize the users' total transmit power 
subject to a set of individual SNR constraints. In these works, the users optimize their uplink power levels and/or 
transmit beams as well as their BS associations. The phenomenon of "cell breathing" is observed whereby the sizes 
of the cells dynamically change according to the congestion levels. Recently references 11251 - 11271 have considered 
the joint BS selection and vector power allocation in OFDM A networks in which the BSs operate on non-overlapping 
spectrum. The users compete with each other for resources in each cell, and at the same time they are able to freely 
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choose to use any (non-interfering) BS in the network. However, it is not clear how these works can be generalized 
to the considered MIMO-HetNet scenario. Reference |[28l is a recent work in such direction. Differently from the 
present work, a downlink HetNet setting is considered. As a result, the user-BS association is determined in a 
centralized fashion, which is very different than the distributed solutions to be presented in this work. 

In this paper, we consider the problem of joint optimization of linear precoder and BS assignment for an uplink 
MIMO heterogeneous network. Our main contributions are summarized as follows. 

1) We establish the NP-hardness of the joint precoder design and user-BS assignment problem for the weighted 
sum rate (WSR) maximization in a MIMO-HetNet. This NP-hardness result is intrinsically different from the 
existing complexity results iTTTl — ifTOTl for the linear precoder design problem where the user-BS assignment 
is fixed. In particular, for the latter results, the NP-hardness of the precoder design problem lies in how to 
determine which users/antennas should be turned off due to strong interference links. However, when the 
user-BS assignment is not fixed, a strong interference link can be turned into a direct link by changing the 
user-BS association, thus effectively mitigating the level of interference. As a result, what makes the joint 
user-BS association and precoder design problem difficult is not the presence of strong interference links, but 
rather its mixed interference pattern. Therefore, the NP-hardness of the joint optimization problem does not 
follow from, nor imply, any of the existing complexity results iPTI- lfTOl for the linear precoder design problem 
when the user-BS assignment is fixed. 

2) We propose a novel game theoretic formulation to find a local WSR-optimal solution of the joint precoder 
design and user-BS assignment problem for a MIMO heterogeneous network. In the proposed game both 
the transmitters and the BSs act as players. We introduce a set of utility functions for the players and show 
that every NE of the resulting game, which gives a set of precoders and user-BS assignment for all users in 
the network, is a stationary solution of the WSR maximization problem. To reach a NE of the game, each 
transmitter greedily determines the best linear precoder as well as the least congested BS for transmission, 
whereas each BS computes a set of optimal prices to charge the users for causing interference. These prices 
serve to coordinate the behavior of the users so that they do not cause excessive interference in the network. 
We show that the resulting distributed algorithm converges to a NE of the game. Notice that the convergence 
of the algorithm implies that the network will be stable in the sense that no user will change BSs indefinitely. 
Simulation results show that this algorithm is very effective in relieving local BS congestion and achieving 
high throughput and load balancing in the network. 

We remark that when the user-BS assignment is fixed, our proposed game reduces to the standard precoder 
optimization game. In this context, the users are again charged with a set of prices that reflect their (negative) 
influence to other users. Thus, our problem formulation can be viewed as a generalization of previous results for 
interference pricing (e.g. [QT|, EES], BHD, llHD) to the MIMO interfering multiple access channel (MIMO-IMAC) 
setting with general system utility functions. 
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The rest of the paper is organized as follows. In Section JIJ we give the system level description of the problem, 
and provide its complexity analysis. In Section IflTl and Section ITVl we formulate the problem into a non-cooperative 
game framework, and study the properties of the resulting games. Furthermore we propose an efficient algorithm 
to compute the equilibrium solutions of the games. In Section |VJ we present numerical results to demonstrate the 
performance of the proposed algorithm. Some concluding remarks are given in Section |VT] 

Notations: For a symmetric matrix X, X y signifies that X is positive semi-definite. We use Tr(X), |X|, 
X^, A(X) and Rank(X) to denote the trace, determinant, Hermitian, spectral radius and the rank of a matrix, 
respectively. We use I n to denote a n x n identity matrix, and use [y, x_j] to denote a vector x with its i th elements 
replaced by y. Moreover, we let M ArxM an d £, NxM denote the set of real and complex N x M matrices, and 
use and to denote the set of N x N Hermitian and Hermitian semi-definite matrices, respectively. Finally, 
we use the notation 0<a_L6>0to indicate a > 0, b > 0, a x b = 0, and use the notation [x] + to represent 
max(x, 0). 

II. System Model, Problem Formulation and Complexity Result 

We consider the uplink of a general MIMO-HetNet consisting of a set M = { 1 , • • • , N} of users that transmit 
to a set Q = {1, ■ ■ ■ , Q} of BSs. Let a be a N x 1 vector representing the system association profile, i.e., a n = q 
means user n connects to BS q. 

Suppose each user n G M has T n transmit antennas and each BS q G Q has R q receive antennas. Let H qjn G 
C R <i xT ™ be the channel matrix from transmitter n to receiver q. Assume T n < R q for all (n, q) G (M x Q), or 
equivalently the channel matrices {H 9 n } are tall matrices. This assumption is reasonable as the number of antennas 
at the BS is typically larger than that of the mobile users. 

Let x n G C T " and y q G <C R<1 denote the transmitted signal of user n and the received signal of BS q, respectively. 
Then y q can be expressed as 

y q = ^ H 9,» x n + Z Q (!) 

neAf 

where z q ~ CN(0, a q Ift q ) is the additive white complex Gaussian noise vector. Suppose there are a maximum of T n 
data streams transmitted by user n and its data symbol vector is denoted by s n G C T ". We assume E[s n s^] = 1^. 
If a linear precoder P n G C T " xT " is used for transmission, then the transmitted vector of user n is x n = P n s n , and 
the corresponding transmit covariance matrix S n is given by S n = E[x n x^] = P n P^ £ §+"• Once the covariance 
matrix S n is computed, the precoder can be obtained by Cholesky factorization. We further assume that each user 
has an individual average power constraint of the form E[Tr(x n x^)] = Tr(S n ) < p n . Define the aggregate transmit 
covariance as S = {S n } n gjv, and the aggregate covariance excluding user n as S_ n = {S m } 

For a fixed association profile a, the interference covariance matrix for user n (at its intended BS a n ) can be 
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expressed as 



C n (S_ n ;a) — a^ n lR an + H ani77l S m H^ 



The covariance matrix of the total received signal at BS q can be expressed as 

N 

G q (S)±a 2 q I + J2H q ,iSiKi- ( 2 ) 

z=i 

For a given user-BS association a, the achievable rate for user n is given by |[30l 

i? n (S n ,S_ n ;a) = log|l Ran +H a „ in S n H^ n C r 7 1 (S_ n )| (3) 
= log I r „ + (sl) if Hf niri C; 1 (S_ n )H a „ )n sl 



log 



(ii) 



loglE-^S)! (4) 



where in (i) we have used the matrix inversion lemma ED; in (ii) we have defined E n (S) as user n's minimum 
mean square error (MMSE) matrix 

E n (S) 4 I Tb - (sl) H H^ >n G a J(S)H a „ in sl >- 0. (5) 

The WSR of the system can be expressed as R(S; a) = YlneM w n R n (S n , S_ n ; a), where the set of nonnegative 
weights {w n }n =1 represent the priorities of different users. We are interested in finding the optimal user-BS 
assignment as well as the transmit covariance matrices that maximize the WSR. This problem can be stated as 

maxi?(S;a) (SYS) 

S,a 

s.t. Tr(S n ) < p n , S n £ S+", VriGjV 
a„ e Q, V n G TV. 

Our first result shows that finding the global optimal solution to the system level problem is intractable in general. 

Theorem 1: Finding the optimal BS association and the transmission covariance matrices (a, S) that solve the 
problem (SYS) is strongly NP-hard. 

Theorem Q] is proved based on a polynomial time transformation from the 3-SAT problem, which is a known 
NP-complete problem [32 j. The 3-SAT problem is described as follows. Given M disjunctive clauses C±, ■ ■ ■ , Cm 
defined on N Boolean variables X± ■ ■ ■ , X^, i.e., C m = t\ V ti V with G {Xi, • • • , X^,Xi, • • • , Xn}, the 
problem is to check whether there exists a truth assignment for the Boolean variables such that all clauses are 
satisfied (i.e., each clause evaluates to 1) simultaneously. We leave the details of the proof to the Appendix lAl We 
note that our complexity result differs from the recent result in |[28l , in which the complexity for the joint user-BS 
association and precoder design problem is analyzed for the downlink direction. 

Motivated by the above complexity result, we focus on designing low complexity algorithms that can provide 
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approximately optimal solutions. In particular, we will consider in the sequel a more general system level problem 
formulated as the following sum utility maximization problem 

max /(S;a) (SYS-U) 

s.t. /(S; a) 4 J2 fn{Rn (S; a)) , R n (S; a) defined in ©, 

neAf 

Tr(S n )<p n , S„eS£", \f n e AT 
a n G Q, Vn£jV 

where / n (-) : — > R is the utility function of user n's data rate. We make the following assumptions on / n (-) : 

A-l) f n {x) is strictly increasing, concave and coercive in x for all x > 0; 

A-2) f n (- log(|X|)) is strictly convex in X for all X y 0. 
Note that this family of utility functions includes well known utilities such as weighted sum rate, proportional 
fairness and the harmonic mean rate utility functions (see lfl"8l ). 

In the sequel, we will develop a distributed algorithm to compute a local stationary solution for the problem 
(SYS-U). Our main approach is based on the noncooperative game theory. 

III. Transmit Covariance Optimization Game for Fixed User-BS Association 

To simplify the presentation of the game theoretic approach, we first consider the case in which the user-BS 
assignment is fixed in advance, and the users are only allowed to optimize their transmit covariances. We will 
design a noncooperative game whose equilibrium solutions correspond to the stationary solutions of the sum utility 
maximization problem (SYS-U). Extension to the general case with flexible association will be presented in the 
next section. 

A. Problem Formulation 

When the user-BS association is fixed, the sum utility maximization problem (SYS-U) can be restated as 

max/(S) (SUM) 
s.t. Tr(S n ) < p n , S n G S*S V n G M. 

Suppose each user fi6jV can optimize its transmit covariance S n . Its feasible set is given by 

Tn = {S„ | Tr(S n ) < p n , S n G • (6) 

Define the joint feasible set of all the users as F = \\ n( zj^F n . 
In order to mitigate interference caused by unintended users, we allow each BS q G Q to post a (matrix valued) 



price Tqn G C Rq to each user n G N. That is, each user n incurs a total charge of Tr J2 q ^Q r ^q,nS n 



for the 
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interference contributed to all the BSs in the network. Define H q = rineAf as ^ e f eas i me set of BS q's pricing 
strategies. Define % = W^qH-q as the joint feasible set of all the BSs. Let T q = {T gn } ng ^, T n = {T q>n } qt =Q, 
and T = {T q } qGQ . 

We model both the users and the BSs as selfish players in a noncooperative game. The players are interested in 
choosing their individual optimal strategies (transmit covariances for the users, and price matrices for the BSs) to 
maximize their own utility functions U n {-) and D q {-). We formulate a covariance optimization game Q c as follows 



A 



{M, Q}, {T,H},{{U n (-)} neAf , {D q (-)} qeQ } 



We need to properly specify the utility functions U n {-) and D q {-) so that the equilibriums of the game Q c correspond 
to the local stationary solution of the sum utility maximization problem (SUM). 

B. The BSs' Utility Maximization Problem 

The BSs' utility functions and their maximization problems are closely related to the structure of the desired 
interference prices. As a result, we start by providing an explicit construction of the interference prices. Let us 
define J\f q = {n : a n = q} as the set of users associated with BS q. Define a m (S) as the derivative of user m's 
utility function w.r.t. to its rate 

df m (x) 



«m(S) 



Ox 



>0 (7) 

x=R m (S) 



where the positivity of a m (S) comes from the condition A-l), i.e., / m (-) is a strictly increasing function. Then at 
a given system covariance S, user n's negative marginal influence to the sum utility of users currently associated 
to BS q is given by 

m<EAf q \n 

= a r?j.(S)H^ Jl C m 1 (S_ rn ) (lf/ q + H gjrn S rn H^ m C m 1 (S_ rn )) H^nSmH^Cj^S-m)!!,^ 

= UfJ J2 o m (S)G- 1 (S)H^„ l s|E- 1 (S)(sl l ) ff H^„G- 1 (S)')H 9 ,, i h (8) 
where in the last equality we have again used the fact that 

1 ~ H H 1 ~ 

E m (S) = I + (Sm) H am Jri C m (S_ m )H amim Sm- (9) 

We propose to mitigate the interference generated by a user n G J\f by charging it a penalty proportionally to 
its negative marginal influence. Specifically, the interference price takes the following form 

T 1,n = ~ Yl V S „/m(^m(S)) h 0. (10) 

m£Af g \n 



s 



This definition of interference price leads to the following definition of a BS g's utility function 

(11) 



D,(T ff ,T_„S) 4 T ,, n - f - £ VsJmiRmiS))) 

n£Af ^ m£Af q \n ' 



Clearly, for a fixed S £ F, the set of prices that maximizes BS q's utility are given by (TTOb . In what follows, 
T 9) n(S) is occasionally used to explicitly indicate the dependency of the prices on the users' transmit covariances. 

C. The Users' Utility Maximization Problem 

To strike a balance between the user's desire to improve its utility and the need to reduce its interference in the 
network, we modify each user's utility as the difference between its true utility and the interference charge 

U n (S n , S- n , T n ) = f n (R n (S n , S_ n )) — Tr [T gjTl S n ] . (12) 

gee 

With this modification, each user's utility maximization problem is given by 

max f7 n (S n ,S_ n ,T n ). (13) 

S„6-F„ 

Note that the function / n (i? n (S n , S_ n )) is a strictly concave function in S n , as it is a composition of a strictly 
increasing and concave function (i.e., / n (-)) and a strictly concave function of S n (i.e., R n (-)) (see IT331 Section 
3.2.4]). As a result, problem ([T3l has a strictly concave objective value and admits a unique solution. In the 
following, we develop an efficient procedure to compute such solution. 

Fix a user n 6 J\f, let q = a n . Let A n = X^/eQ^i," — 0- Using these notations, user n's utility maximization 
problem ([T3l can be written as 

max f n (log |I + H^SnH^C- 1 !) - Tr[A n S n ]. (14) 

The Lagrangian of this problem is given by 

L(S n , fi n ) = f n (log |I + H^SnH^C- 1 !) - Tr[(A n + Mn I)S n ] + ^ n p n (15) 

where fi n > is the Lagrangian multiplier for the power constraint. The dual function is d(/j, n ) = maxs„^o L(S n , fj, n ) 
The optimal primal-dual pair (S*,/x*) should satisfy the KKT optimality conditions 

s n = arg max L(S n ,/i*), S* ^ (16) 
0<<±(p n -Tr[S;])>0. (17) 

For any fixed fj, n > 0, the solution to the problem argmaxg„ej-„ L(S n ,fJ, n ) can be obtained as follows. Using 
the fact that A n ^ 0, then for any fj, n > 0, we can perform the Cholesky decomposition A n + fj, n I = L^L n , which 
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results in Tr[(A n + /i n I)S n ] = Tr[L n S n L^]. Define S n (/x n ) = L n S n L^, then we have 

L(S n ,fj, n ) = f n (log |C„ + H qjn L~ 1 S n (iJ ln )L~ II H^ in \ + c x ) - Tr[S„(/i„)] + fi n p n 

( = } /„ (log |I + B.- 1 H,,„L,; 1 S n ( M „)L r ; ff Hf n B r ; H | + c 2 ) - Tr[S„( M „)] + /i„p„ 
( = } /„ (log |I + F n A„M r f S„( M „)M„A n F^| + c 2 ) - Tr[S„( M „)] + fi n p n 
= fn flog I + A„S„(^„)A„ + c 3 ) - Tr[S„(/i n )] + ^„p„ = L(S„(/i n )) (18) 



where ci, C2, C3 are some constants that are not related to S n . In step (i) we have used the Cholesky decomposition: 
C n = B n B^; in (ii) we have used the singular value decomposition B^Hq^L" 1 = F n A n M„ ; in (iii) we have 
defined S n (/%) = M^S n (/i n )M n and used the fact that M n and F n are unitary matrices. Note that if fi n = 
and A n is not full rank, we can use generalized inverse to replace L" 1 . We then argue that the (unique) optimal 
solution S*(/i n ) to the following problem must be diagonal 

max L(S n (p n )). (19) 

S„(/i n )X0 

Assume the contrary. Note that we have I + A n S n (/U n )A n y 0. Then from the Hadamard inequality ||3T| , we can 
always remove the off-diagonal elements of the optimal solution and increase the value of log I + A n S* (/i n )A n 
while keeping the value of Tr[S* (/i n )] unchanged. Since f n (x) is a strictly increasing function, the objective value 
is also increased, a contradiction to the optimality of the non-diagonal solution. 

Let Si = [S n (nn)]i,i- Utilizing this diagonality property, solving the matrix optimization problem (fT9l reduces to 
a vector optimization problem of the form 

(T n \ T n 

log(l + Si[A n ]?J + c 3 )-Y j s i (20) 
i=i / i=i 

Let denote the Lagrangian multiplier associated with the constraint Sj > 0. The optimal primal-dual variables 

{s*,Q}J=i must satisfy the following optimality conditions 



dfn(x) 



dx 



?. 

1,1 



2 



V » = !,-•• ,T n (21) 



0<Ci<>0, Vi = l,--- ,T„. (22) 
The condition ((211) implies that there must exist a constant c* that satisfies 

c *4 l + i4) (1 - c) ' Vi ' 

where c* = «J Ef=il°g (l + <[An]-,i)j, with a n (i? n ) = f^M^^. Due to we have 

c * = Q "(l log (^§f))' 
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Suppose we have the optimal c* (note that due to the strict concavity of problem (|20l i, c* is unique), then using 
condition (1221 and the definition of c* in (|23T ). the optimal primal-dual variables {s*,C*}^i can be expressed as 

C ~ IA^~' llc - |A^~ > U 



c 
0. 



11 c ~ ia^: 

otherwise, 



0, 



if c* 



;i-c*[A n ] 



2 ) 



otherwise. 



> 



(24) 



In short, s* = [c* 



andC= [l-c*[A4 



Our task then becomes finding the optimal c*. For any c > 0, let us define (*(c) = [l — c[A n ]?^ . We have the 
following lemma characterizing the relationship between a n (-) and c. The proof of it can be found in Appendix 151 



Lemma 1: For any fixed c, if a n ( Yli=i^°& (t^T^Tct) ) ^ c ' ^ nen c * > c - Otherwise, c* < c. 



Q(c)> 

This lemma suggests that we can find the optimal c* by a bisection search. 

TABLE I 

The Binary Search Procedure to Find c* 



51) Choose c u and c l such that c* lies in [c l , c u j. 

52) Let c mid = (c ; + c M )/2 



Compute s*(c mid ) 



[c mid 



[An]? t - 



and C(d 



mid\ 



[1 



„mid 



Mi 



S3) If aJ E£ilog(Y 



„mid 



-C*(c mid ) 

S4) Go to S2) until the desired accuracy is reached. 



> c mid , let c l = c mid . Otherwise, let c 11 = c 



mid 



Note that in the special case where the utility function is the weighted sum rate utility, i.e.,f n (R n ) = w n R n , we 



V i = 1 



,T„. 



have c* = w n . Hence no bisection search is needed, and we directly obtain s* = 
This is the well-known water-filling solution. 

Once we have the solution S*(/t n ), we can obtain S*(/i„) = L~ 1 M n S* (/i„)M^L~ H . Note that L(S n ,/i n ) is 
a strictly concave function of S n . As a result, for a fixed \i n , the solution S* (/i n ) is unique. 

The remaining task is to find the optimal /x* that satisfies the complementarity condition (flTT i. From a general 
result on penalty method for optimization, e.g., ||34l Section 12.1, Lemma 1], the solution Tr[S*(/i n )] must be 
monotonically decreasing with respect to [i n . Such monotonicity result suggests that we can find the optimal /x* 
that satisfies the complementarity and feasibility condition (ITTb by the following bisection search procedure. 

TABLE II 

The User' s Utility Maximization Procedure 



51) Choose fj% and fi l n such that fi* n lies in [fj, l n , 

52) Let /x™ d = (fi l n + Hn)/2. Compute decomposition: 
A n + /J. mid I = L,^L n 
C n = B n B n 

B^Hg^L" 1 = F n A„M^. 

53) Compute c* using the procedure in Table U 

54) Compute S;(/C d ) by [S*(^ id )] M = [c* - + , 

55) Compute S;« ld ) = L^M^/C^M^L-^.'* 

56) If Tr(S;(/i™ d )) > p„, let = ^ id ; otherwise let 

57) If |Tr(S;(/i™ d )) -p n \ < e or - fj, l n \ < e, stop; otherwise go to S2). 



1- • ,T n . 



, .mid 
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D. Analysis of Nash Equilibriums (NE) 

Consider the game Q c . We first show that our choices of the utility functions give rise to a nice relationship 
between the utility of the users and the system utility function f(S). 

Proposition 1: If the conditions A-l) and A-2) are satisfied, then we have 



(S))-U n (S (S)) > 

^ /(Sri) S_ n ) — f(S n , S_ n ) > 0, V S n , S n G J 7 n, S„ n G J~—n- 



(25) 



Proof: Fix any S € J 7 and S € J 7 . Pick an m / n. Utilizing the assumption that / m (— log |X|) is convex in 
X, we can linearize f m (— log |E m (S n , S_ n )|) at the point E m (S n , S_ n ) using Taylor expansion 



/ m (-log|E m (S„,S-n)|) > / m (-log|E m (S)|)+Tr V Em / m (E m )(E m (S n ,S_ n ) - E m (S)) 



/ m (-log|E m (S)|) 



dfm{x) 



d.v 



Tr 



E- 1 (S)(E m (S„,S_„)-E m (S)) 



(26) 



a=-log|E m (S)| 

To proceed, we need the following lemma whose proof is relegated to Appendix ICl 
Lemma 2: The function —Tr E r " n 1 (S)E m (S n , S_„) is convex in S n , for all n ^ m. 

This convexity result allows us to perform a further linearization in the variable S n via Taylor expansion of 
E m (S n ,S_ n ) at S n 



- Tr 
> -Tr 
= -Tr 



E" 1 (S)(E m (S 

-^m 1 (S)(Sm)^H^ ifn G fn 1 (S)H ai7iin (S n — S n )H^ njn G fn 1 (S)H ai7i/m Sf 
H a / „,nG m 1 (S)H amim SmE m 1 (S)(Sm) H H a ^ im G m 1 (S)H ami „(S n — S n ) . 

Plugging (l27l) into (l26l) . utilizing the expression for the negative marginal influence in (|8), we obtain 

J2 /m(-Rm(S n ,S_ n )) = J2 /m(-l0g|E TO (S„,S_„)|) 

> |/ m (i?m(S))+Tr[v s „/ m (i? m (S))(S„-S ? 



(27) 



This inequality implies 

N 

/(S n ,S_ B )-/(S)=/ n (S 

m^n m— 1 

> / n (i? n (S n ,S_ n )) + ^ [Tr [v Sn / m (Ji ro (S))(S„ - §„)] ] - /„(i?„(S)) 

= fn{Rn{Sn, S_„)) - J! Tr [T g , n S„] - I f n (Rn$)) ~ X) Tr [T 9 ,„S„] 
<?es \ <?eQ / 



(S)) - J7„(S (S)). 



(28) 



This inequality implies that (1251 ) is true. 
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The property (|25l l is essentially the generalized potential property Q for a class of so called Potential Games |[35l , 
with only one subtle difference that in (T25T ) the implication is dependent on a "state variable" T n (S). 

Using Proposition [T] we can establish the one-to-one correspondence between the pure NE points of the game 
Q c and the KKT points of the sum-utility maximization problem (SUM). Recall that a pure strategy NE of the 
game Q is a tuple of strategies {S*,T*} such that the following set of inequalities are satisfied 

U n (S* n ,S*_ n ,T* n ) > U n (S n ,S*_ n ,T* n ), V S n G T n , V n e M 
D n (T;,T!L,,S*) >£> n (T„T*,,S*), VT,£?i 9 ,VgeQ. 

By utilizing Proposition [T] we have the following characterization of the NEs of the game Q c . 

Theorem 2: The tuple (S*,T*) is a NE of the game Q c if and only if S* is a KKT point of the problem (SUM). 
Proof: We give an outline of the proof here. First suppose (S*,T*) is a NE of the game. Then from the 
definition of NE we have that for any S n G F n and n G TV, £/ n (S* , S* n , T* ) > f7 n (S n , S* n , T* ). Using (l25i 
we have R{S* n , S*J > R(S n , S*J for all S n G F n . This means S* = argmax s „ e j„ R(S n , S*_ n ), V n G TV. We 
can verify that the KKT condition of this N problems is the same as the KKT condition of the original (SUM) 
problem. The other direction can be proved similarly. ■ 

IV. Joint BS Selection and Transmit Covariance Optimization Game 

In this section we extend the game theoretic framework described in Section [III] to the case where the user-BS 
associations are not fixed. 

Let us define user n's joint strategy as J n = (S„,a n ), and its feasible space as J n = F n x Q. Let J_ n = 
(S_ n ,a_„), and J = {J n }n£j\f- In this case, each user's rate is still defined by ©, but we have to make the 
dependency of association profile explicit. We use R n (J n , J_ n ) to denote user n's rate. We use C n (S_ n , [g, a_ n ]) 
to denote the interference covariance that user n would have experienced if it selects BS q, while all other users 
use the strategy (S_ n ,a_ n ). Let N q (&) denote the set of users associated with BS q under association profile a. 
Moreover, to make the dependence of the sum utility maximization problem on the underlying user-BS association 
explicit, we use SUM(a) to denote the sum utility maximization problem when the association profile is chosen as 
a. 

Let U n (-) an d Aj(0 denote user n and BS q\ utility functions, respectively. The joint BS selection and covariance 
optimization game Q J is defined as 

G J = | {TV, Q},{j,n},{{U n (3,T n )} neU ,{D q (T,3)} qeQ }y 

'The generalized potential property is referred to as the following relationship between the players' utility functions and a "potential 
function" P(-): let x n be player n's action profile; for any two x n ,x„ G Xn< f° r all x_ n G X-n, and for all player n, [7n(x n ,x_ n ) — 
Un (x^ , X_ n ) > implies P(x n ,x_„) — P(x ) > 0. 



13 



We refer to the game Q J as a hybrid game, because the strategies of a subset of the players consist of a covariance 
matrix and a discrete index. We define the utility functions £7 n (") and D q {-) similarly as in ([T2l and (fTTb 

U n (J n , J_ n , T n ) = / n (i? n (J n , J_ n )) — ^ ' Tr [Tg ;W S n ] 
Aj(T g ,T_ 9 , J) = T 9 , n - (~Y1 Vs./ m (H m ((S B ,a B ),J_ m ))j • 

nG-V m6A/',(a)\ti ' 

Note that both of the utility functions defined above are dependent on the user-BS association vector a. In order to 
emphasize the relationship between the optimal solution of BS q and the users' strategies J, we occasionally use 
Tg,n(J) or T 9 n (S,a) (resp. T(J) or T(S,a)) to denote the optimal prices charged by BS q to user n (resp. the 
set of prices charged by all the BSs). 

The pure NE of the game Q J is the tuple (J*, T*) that satisfies 

u n (r n , r n T*n) > u n (Jn, J* n ,T*) , v 3 n e j n , v n e M 

D q (T* q ,T%,J*) > D q (T q ,T*_ q ,3*), V T q e^.VgeQ. 

The game Q J with the utility functions defined above again possesses the "generalized potential" property, which 
is essential in establishing the correspondence between the pure NEs of game Q J and the stationary solutions of 
the sum-utility maximization problem. 

Proposition 2: For any utility function f n (') that satisfies the assumptions A-l) and A-2), we have 

U n (j n , J_„T n (J)) - U n (j„, J_„T,„(J)) > (29) 
^* / ( Jfn J— n ) / ( Jfi) J— n ) > 0, V J n , J n G J7ni J — n G 



The proof of this proposition is relegated to Appendix |D] The key observation used in the proof is that 

-Rm((Sn> a n ), J _ n ) = R m {(S n , a n ), J _ n ), V rn ^ n 
T g , n (S,a) = T ? , n (S, [a„,a_ n ]), V<jGQ. 

That is, if user n unilaterally switches from BS a n to a n but keeps its covariance matrix unchanged, then all other 
users' transmission rates as well as the price charged for user n remains the same. 

Due to the hybrid structure of the users' strategy space J, conventional existence results of the NEs for a 
iV-person concave game (e.g., (36]) do not apply here. Fortunately, by utilizing Proposition |2j we can extend our 
argument in the proof of Theorem |2] to show the following existence result of the NE of game Q J . 

Theorem 3: The game Q J must admit at least one pure NE. Moreover, if (S*, a*, T*) is a NE of the game Q J , 
then S* must be a KKT solution of the problem SUM (a*). 

Proof: We first claim that the global optimal solution of problem (SYS), say (S,a), along with the corre- 
sponding price matrices T(S,a) is a NE of the game Q J . Assume the contrary, then there must exist a user n € N 
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with (S n ,a n ) ^ (S n ,a n ) who has incentive to switch 



Un y(S n , a n ) , J—n, X n (J) J > U n ^J n , J_ n ,X n (J 
However, using Proposition |2] this implies that 



(30) 



(31) 



which contradicts the global optimality of (S, a). The second part of the theorem can be shown following the same 
proof in Theorem |2] ■ 
In the following, we propose a distributed algorithm for the users to reach a NE of the game Q J . Central 
to the proposed algorithm is the procedure developed in the previous section that allows the users to compute 
their transmit covariance matrices. The algorithm works by alternating between the users' and the BSs' utility 
maximization problems. In each iteration t, a single user n G M updates its transmit covariance and BS association 
by solving its utility maximization problem 



(S^X) = arg max U n ((S n , a n ), Jl" 1 , T^" 1 ) . 

(S„,a„)ej"„ 

Then all the BSs update their interference prices by solving their respective utility maximization problems 



(32) 



arg maxDg(T g , T 



t-i 



VgGQ. 



(33) 



Each user's utility optimization problem (1321 1 can be performed by the following two steps: a) solve Q inner 
covariance optimization problems maxs„ e j„ U n ((S n , q), J_ n , T n ), one for each BS q G Q (each of these problems 
can be solved using the procedure listed in Table HT1 in Section UTTl): b) pick the best BS in terms of the optimal 
value of the inner covariance optimization problem. 

The proposed algorithm naturally incorporates the joint optimization of BS association and linear precoder into 

individual mobile users' optimization problem. The detailed algorithm is listed in Table UTTl 

TABLE III 

A Distributed Best-Response Algorithm 



51) Initialization: Let t = 0, each user n G M randomly choose G J n \ 
each BS q G Q set T° = 0. 

52) Choose a user n G M. Compute J2 q eQ Tq,n- 

53) User n computes J^ +1 by solving: 

(Sn +1 X +1 ) = argmax 96Q max s „ 6 jr i U n ((S n , q), Si n , a*_ n T^). 
For the rest of users m ^ n, set J^ 1 = 3 l m . 

54) Each BS q G Q updates its price matrices by 

Tqtn = ~ Sm.eA/' g (a t + 1 )\n V S„ fm(R m (S„ 1 , S_ m )), V 71 G M. 

55) Continue: Set t = t + 1, go to S2) unless some stopping criteria is met. 



An important feature of the algorithm is that the computation of each of its steps is closed form (subject to 
efficient bisection search) and distributed. We briefly remark on the distributed implementation of the algorithm. 
The following three assumptions are needed for this purpose. First, local channel information is known by each 



15 



user, that is, each user n has the knowledge of {H gjn } g6 g. Second, each BS has a feedback channel to all the 
potential users. Third, each BS knows the utility function of the users that are associated with it. 

Under these assumptions, the proposed algorithm can be implemented distributedly. At iteration t, each BS q can 
compute the price T* locally by measuring the total received signal variance G q (S t ) and computing the MMSE 
matrix E m (S*) of each associated user m £ M q (cf © and ©). Suppose user n is scheduled to update at iteration 
t. Then all the BSs can broadcast their pricing information T* for user n as well as their total received signal 
covariance G 9 (S*) (note, due to symmetry, only upper triangular parts of these matrices need to be transmitted). 
Upon receiving all this information, user n can carry out its utility maximization locally. 

Moreover, when the network is operated in a time division duplex (TDD) mode, the information that needs to 
be broadcast can be significantly reduced. This reduction is made possible by utilizing the following two facts: 
1) in TDD mode, the uplink channels can be viewed as the Hermitian transpose of its reverse channels (i.e., 
Hq,n = H^ g ); 2) each user n only needs the sum of all the prices charged for it: XLeQ ^-q,n- Specifically, the BSs 
do not need to broadcast the pricing information for user n explicitly. Each BS q only needs to broadcast G ? (S*) 
by using the following transmit covariance matrix 

E « m (S)G~ 1 (S)H g , m S,iE m 1 (S)(si) H H^ m G- 1 (S). (34) 

m£Af,\n 

This matrix can be calculated once BS q obtains the measurement of the total received signal variance G ? (S*) 
and computes the MMSE matrix E m (S*), for all m G J\f q (a t ). In this way, user n can decode the messages and 
measure the total received signal covariance expressed as 

Y,Kn\ E " m (S)G- 1 (S)H,, m S^E- 1 (S)(si) if H^ m G- 1 (S)) H,, n . (35) 

qeQ \m£Af q \n J 

By the definition of prices in (H) and (fTOb . the received signal covariance is precisely the total price YlqeQ^-q^- 
In the following we provide the convergence result of the proposed algorithm. The details of the proof are 
presented in Appendix |E] 

Theorem 4: The sequence {/(S*, a')}^ generated by the proposed algorithm is monotonically increasing and 
always converges. Any limit point of {S*, a*, T*}^ is a NE of the game Q J . 

We remark that the proposed algorithm can also be applied to the scenario in which the user-BS assignment is 
fixed. In this case the users only need to perform a single inner optimization in S3). An immediate consequence 
of Theorem |4] is that this reduced form of the algorithm also converges to the NE of the game Q c . 

Corollary 1: When the user-BS assignment is fixed, the sequence {/(S generated by the proposed algorithm 

is monotonically increasing and converges. Any limit point of {S*, T*}^ is a NE of the game Q c . 

This corollary generalizes the convergence result presented in (371, which deals with only the single antenna 
case and is limited to the weighted sum rate utility. Furthermore, it establishes the convergence for the algorithm 
proposed in ifTTI . 11121 . since the latter is a specialization of our algorithm to the case where the user-BS association 
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is fixed and the system utility is the weighted sum rate function. 

V. Numerical Results 

In this section, we compare the performance of the proposed algorithm with the WMMSE algorithm ||18:j. The 
latter is known to be an effective method to solve the sum utility maximization problem for the MIMO interference 
channel, except that it requires the user-BS assignment to be fixed. To facilitate the comparison, we fix the user-BS 
assignment for the WMMSE algorithm using the received signal strengths, as is done in the conventional cellular 
networks. We demonstrate that the distributed algorithm proposed in this paper can achieve a higher spectrum 
efficiency and more effective load balancing in a HetNet than the WMMSE algorithm. 

We consider a single macro cell in a HetNet. The macro cell consists of 7 pico cells, each containing 1 pico BS, and 
has a total of 16 users. The distance between adjacent pico BSs is 200 meters (representing a dense macro cell with 
small pico cell sizes). Let d q:1l denote the distance between pico BS q and user n. The entries of the channel H g n 
are generated from distribution £/V(0, a* n ), where the standard deviation is given by a q ^ n = (200/d gjn ) 3 ' 5 L q , n , 
and 10 log 10(L gjn ) ~ AA(0, 8) is used to model the shadowing effect. We fix the environment noise power as 
o q = 1,V g £ Q, and let all users have the same transmit power limit p n = p, V n G M. We define the signal to 
noise ratio as SNR = 10 log 10 p. 




Fig. 1. Comparison of the performance of the proposed algorithm and the WMMSE algorithm when the proportional fair utility is used. 
N — 16, Q — 7, users are all located at the cell edges. Left: comparison of the users' averaged rates. Right: comparison of the CDF of 
users' achieved rates. 

We first consider a scenario in which the users are all located at the pico cell edges, and one of the pico 
BSs is congested. In particular, we place half of the users uniformly at the cell edges of BS 1, which is within 
di >n € [90, 100] meters. We place the rest of the users randomly at the cell edges of other pico BSs. For the 
WMMSE algorithm, we let the users associate to the pico BSs with the strongest direct channel (in terms of the 2- 
norm of the channel matrices). For our proposed algorithm, we place a restriction that the users can only choose their 
association among the three strongest pico BSs. We initiate our algorithm by assigning the users to their respective 
strongest pico BSs. Fig. Q] compares the performance of the algorithms when the proportional fairness utility is 
used, i.e., f n {Rn) = log(-Rn), V n G M. Each point on this figure is averaged over 100 randomly generated user 
positions and channel coefficients. The left panel of Fig. Q] compares the users' averaged rates achieved by different 
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algorithms. The right panel of Fig. [TJcompares the CDF (cumulative distribution function) of the individual rates of 
the two algorithms when SNR = 30dB. Fig. Q] shows that if the user-BS assignment is allowed to be optimized, the 
proposed algorithm can achieve a substantially higher spectrum efficiency and fairer rate allocation, as compared to 
the case when the user-BS assignment is fixed. This is reasonable since assigning weak users to less congested BSs 
(rather than the closest BSs) effectively reduces the interference level (hence the congestion level) of the congested 
BS. In this way, both user fairness and the spectral efficiency of the entire network are enhanced. 




SNR(dB) Individual Rate 



Fig. 2. Comparison of the performance of the proposed algorithm and the WMMSE algorithm in a HetNet with different BS capabilities. 
The proportional fair utility is used. N = 16, Q = 7. Half of the users are uniformly distributed in cell 1, all other users are uniformly 
distributed in other cells. Three neighboring BSs have 10 receiving antennas, and all other BSs have the same number of antennas as BS 1 
(which is either 2 or 4). Left: comparison of the users' averaged rates. Right: comparison of the CDF of users' achieved rates. 

To highlight the load balancing capability of the proposed algorithm, we next consider a scenario in a HetNet 
where the pico BSs have different capabilities. Specifically, three out of the six neighboring pico BSs of BS 1 
have 10 receive antennas, and all other pico BSs (including BS 1) have fewer receive antennas. Half of the users 
are uniformly located in cell 1 (which is within d\ >n G [20, 100] meters) and the rest of the users are uniformly 
located in other cells (d q>n G [20, 100] meters, q / 1). Again we use the proportional fairness utility function, and 
the proposed algorithm compares favorably with the WMMSE algorithm. See Fig. [2] 

Interestingly, by using the proposed joint BS association and linear precoder optimization algorithm, the "cell 
breathing" phenomenon GUI can be observed. This phenomenon refers to the desirable load balancing property of 
a network: when a cell is congested, it contracts and the cell edge users automatically switch to adjacent cells. See 
Fig. [3] for an illustration. 

VI. Conclusion 

In this paper, we consider the joint design of the user-BS assignment and the users' linear precoder in a 
multicell heterogeneous network. By a careful user-BS association, users in a hot spot can avoid congesting the 
nearest BS or causing excessive interference to each other. Unfortunately the overall joint optimization problem is 
shown to be computationally intractable. To find a high quality locally optimal solution, we propose an efficient 
and low-complexity algorithm using a game theoretic formulation. The effectiveness of the proposed algorithm 



18 





Fig. 3. Illustration of the cell breathing phenomenon. Q = 7, N = 20, R n = 4, T q = 4, proportional fair utility is used. In the figures, the 
diamonds denote the BSs, the dots denote the users, and the lines indicate associations. Left: user-BS assignment in which users are associated 
to the BS with the strongest channel magnitude in terms of 2-norm. Right: user-BS assignment generated by the proposed algorithm. 

is demonstrated via numerical simulations which show substantially improved spectrum efficiency and fairness 

provisioning. A drawback of this algorithm is the fact that it requires the exact knowledge of channel state 

information (CSI). An important issue worth investigating is to what extent we can use inexact channel state 

information or just use the long-term channel statistics in place of the CSI. Another interesting issue is to incorporate 

the users' quality of service constraints in the problem formulation. 

Appendix A 
Proof of Theorem Q] 

Proof: Given M disjunctive clauses C±, ■ • ■ , Cm defined on N Boolean variables X\ ■ ■ ■ , Xn, i-e-, C m = t\ V 
*2 Vt3 with ti £ {^i, • • • , Xjy,Xi, • • • , Xn}, the 3-SAT problem is to check whether there exists a truth assignment 
for the Boolean variables such that all clauses are satisfied (i.e., each clause evaluates to 1) simultaneously. Let 
7r(C m ,i) denote the i th term of the clause C m , and let X(t) denote the index of a term tfs corresponding variable. 
For example if C m = XiV X 2 V X 4 , then 7r(C m , 1) = X u and Z(7r(C m , 3)) = 4. 

Given any 3-SAT problem with M disjunctive clauses and N variables, we construct an instance of multiple BS 
multi-user uplink network with 3M + TV BSs and M + 2N users. Let T n = R q = 1, cr^ = 1, p n = 1 and w n = 1, 
for all q, n. Let h q , n denote the channel coefficient between user n and BS q. For each clause C m , we construct 3 
clause BSs {c^, c^, c^}, and construct 1 clause user denoted as user C m . For each variable X n , we construct 1 
variable BS denoted as x n , and construct 2 variable users denoted as X n , X n . The channel coefficients are specified 
as follows. The clause user Cfti hcts nonzero channels only to the clause BSs {c l m }f =1 . The variable user X n has 
nonzero channels only to the variable BS x n and the clause BS C* m that satisfies X n = 7r(C m , i). Similarly, the user 
X n has nonzero channels to the variable BS x n and the clause BS C % m that satisfies X n = ir(C m ,i). Specifically, 
the channel coefficients are designed as: 



h, 



q,C v 



a/7 a G ic 1 c 2 c 3 ) 



(36) 
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l q,X n 



V?, q = %n, 

1, q = d ml and 7r(C m , i) = X n , 
0, otherwise. 



(37) 



\/7, q = x n , 
h q ,X n = S 1, g = C and 7r(C m ,i) = X n , 
0, otherwise. 

To illustrate, for a given clause C m = X\ V X2 V X3, we construct the network shown in Fig. @] 



(38) 




o 

*i x 2 x 2 x . 



.3 .1.3 

c q = x,vx 2 vx 3 \ = 4l h 2 = 1 /z, = V'7 

Fig. 4. Construction of the network for clause C m 




1 1.5 2 

Achievable rate of user 1 



XiVX 2 VX 3 . 



Fig. 5. Achievable rate frontier of two-user single channel 
interference channel when all the direct and cross talk channel 
coefficients are set to be v7. 

Our claim is that the 3-SAT problem will be satisfied if and only if the network we constructed achieves a total 
sum rate of at least 3(M + N). 

Suppose that the 3-SAT problem is satisfied, then we perform the following assignment: 1) If X n = 1, assign 
the corresponding variable user X n to BS x n . Otherwise, assign user X n to BS x n . 2) for each clause C m , pick 
an index i* m G {1, 2, 3} such that ir{C m ,i* m ) = 1 (note that because the 3-SAT problem is satisfied, we can always 
do so). Assign the clause user C m to the clause BS Cm . We claim that by the above user-BS assignment, and 
by letting all the assigned users transmit with full power, the overall sum rate achieved is 3(M + N). To argue 
this claim, we first note that there is a single user in the variable user pair X n , X n that transmits with full power 
to BS x n . Thus each variable BS x n is free of multiuser interference and obtains a rate of log(l + \) = 3. We 
then consider an arbitrary clause C m , and pick a term t that evaluates to 1, i.e., t = -K{C m ,i* m ). According to our 
assignment scheme, user i does not transmit while user t transmits with full power. By our construction of the 
channel in (I36l)-(l38l). the only variable user that has nonzero channel to the clause BS Cm is user i. Since user 
i does not transmit, then the clause user C m can transmit to clause BS 4' free of interference. Consequently it 
obtains a rate of log(l + |) = 3. In summary, each variable BS is able to achieve a rate of 3, while each set of 
three clause BSs {c^, c^, c^} obtains a total rate of 3. Thus the total system sum rate is 3(M + N). 

Conversely, suppose the network achieves a total rate of 3(M + N), we argue that the corresponding 3-SAT 
problem must be satisfied. We show this direction by three steps. 
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Step 1) We first claim that for any variable BS x n , its maximum achievable sum rate is obtained when a single 
variable user, X n or X n , transmits to it using full power, while the remaining user does not transmit. Let px n and 
Px denote the transmission power of user X n and X n , respectively. The rate region of this uplink channel when 
the interference is treated as noise can be expressed as 



TZ = {(R Xn ( PXn ,PxJ,RxAPx n ,PxJ) -o <Px n <1,0<P X „ <i} 

l x n ,X„ 



where R Xn (px n , P Xn ) = lo S 1 + 



PxJha 



1 +PxJK n ,x, 



2 



Rx n (Px n ,PxJ = log ( 1 + P ^ K :^ - ) (39) 

with the channel coefficients given as h x x = hx n ,x n = \/7. In IT381 . Charafeddine and Paulraj derived a complete 
characterization of this rate region (without time-sharing operation). This region (plotted in Fig. |5]l asserts that the 
sum rate maximum point can be achieved only if a single user transmits using its full power (at point A or B). 

Step 2) We then argue that a variable user should never be assigned to any clause BS at the sum-rate optimal 
solution. To see this, consider an arbitrary C m and an arbitrary i G {1,2,3}, let t = n(C m ,i). Suppose at the 
optimal solution user t transmits to BS c\ n , then it can obtain a maximum rate of log(l + j) = 1. If BS xx{t) nas 
no associated user at the optimal solution, then the same user t can switch to BS xx(t) an d obtain a rate increase 
of log(l + f) - log(l + 1) = 2. This is a contradiction to the optimality of the solution. Consequently BS x X {t) 
must have user i associated to it. In this case, the maximum sum rate that user t and i can obtain is the optimal 
solution for the following problem 



max log ( 1 + 7Pt ) + log (1 + p t ) . 



0<Pt,Pt<l 

Clearly at the optimal pi = 1, and the tuple (pt,Pt) = (0,1) is a feasible solution to the above problem with 
an objective value log(l + V) = 3. We then argue that (pt,Pt) = (0,1) is in fact the optimal solution to this 
optimization problem. Specifically, we will show that the following is true 

log (l + —L— ) + log (l+p t ) < log(l +7), V < p t < 1. (40) 



l + 7 Pt 

Equivalent we show the following inequality 



+ TTWt ) (1 +Pt) < 8 ' V ° < Pt ~ L (41) 

Note that = 1 — (x+7p M • So when p t G [0, the slope of f(pt) is negative, and the functional value of 

f(pt) is strictly decreasing, when p t G ( v ^~ 1 , 1], the functional value of f{pt) is strictly increasing. Combining 
with the fact that /(0) = 8 and f(l) = 2 + \ < 8, we have that f(p t ) < 8 for all < p Xn < 1. In summary, a 
variable user never transmits to any clause BS at the sum rate optimal solution. 

Step 3) The third step is to show that any clause BS tuple {c^, c^, c^} can at most obtain a rate of log(l+ j) = 3 
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at the sum rate optimal point. This step is trivial as the only candidate that can select and transmit to clause BSs 
{c^j}f =1 is the clause user C m . When there is no interference, the maximum rate user C m (and equivalently the 
set of BSs {cj 7l }? =1 ) can get is log(l + 7) = 3. Note that this maximum rate can be achieved only when at least 
one of these clause BSs does not experience interference from the variable users. 

Step 2-3 together imply that each variable BS and each clause BS tuple {c^}? =1 are able to achieve a rate 
at most 3 at sum rate optimal point. As a result, in order to achieve a system rate of 3(M + N), each of them 
must achieve a rate that exactly equals to 3. From Step 3, for any clause C m , if a BS tuple {c^}?^ achieves a 
throughput of 3, then it is only possible that there exists a index i* m and a user t* m = C(m,i* m ) such that user t* m 
does not transmit. Set the terms t* m = 1 for each clause m = 1, • • • , M. This assignment ensures that every clause 
C m contains at least one term that is assigned to 1. Consequently, the 3-SAT problem is satisfied. 

Since the above construction only involves universal constants, the joint linear precoder design and base station 
selection problem (SYS) is strongly NP-hard. The proof is complete. ■ 

Appendix B 
Proof of Lemma Q] 

Proof: We prove this claim by contradiction. Suppose the contrary, that a n ^Y^i=i^°S ( (' c ) ) 1 > c anc ^ 
c* < c. Define s*(c) = [c — nA? ] + . Define the following two sets 

l(c)±{i\s*(c) = 0}, X(c)4{i| a *( c )>0}. 
From the definition s|(c) = [c- , A \i ] + , C*( c ) = i 1 ~ c [ A n]i,i] + > we have 



\iex( c ) 

Note that if c* < c, we have X(c*) C 1(c), X(c*) 5 1(c). This together with the fact that c* < c imply 

53 log(c[A re ]y > |3 log(c*[A n ]^). 

iex(c) «ex(c*) 

Since, as the derivative of a concave utility function, the function a n (-) is monotonically decreasing, we obtain 



v iei(c*) 



>a„ ( 53 log(c[A„]y j >c>c*. 

iGX(c) 
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This is a contradiction to the optimality condition that 



a, 



c . 



Thus we conclude that any fixed c, if a n ( J2i=i l°g {t -c(c) , 
The other direction can be shown similarly. 



■c\A r , 



> c, then c* > c. 



Appendix C 
Proof of Lemma [2] 



We show that the function — Tr 



E m 1 (S)E m (S n , S_ n ) 



is convex in S„. Let 



5m(S n ) — Tr E m i (S)(Sm) i ^H^ G fn i (S)H ami?Tl S 



From ll33l . we see that in order to prove that <7 m (S n ) is convex in S n , it is sufficient to prove that g m (S n + tD) 
is convex in the scalar variable t, for any fixed direction D G S Tn as long as S n + iD G ST™. In what follows, we 
will show that for all D G § T ™ that satisfies S„ + tU y 0, we have gM|i+*g) > o. 

To proceed, we make the following definitions 

B — H a?i „DH^ n , 

JV 



The first and second order derivatives of <? m (S n + iD) with respect to t can be expressed as 
dg m {S n +tT>) 



dt 

d 2 g m (S n + tT>) 



dt 2 



-Tr 



= 2Tr 



E- 1 (S)(S= l ) ff Hf mim G- 1 WBG- 1 WH am , m S^ 



E t n 1 ( s ) ( s "> ) H H f„ , m G,,/ (t)BG m 1 (t )BG Tn 1 (i )H am ,mS, : 

The fact that S„ + iD y ensures G m (t) y 0, which further implies 

(Sn)^H^ n G m 1 (t)BG~ 1 (i)BG~ 1 (t)H arein S^ ^ 0. 

Combining with the fact that E^S) >- 0, we conclude that ^(s^+m) > . 

Appendix D 
Proof of Proposition [2] 

Proof: We first write an equivalent form of (|29) (note that we have defined J = (S,a)) 



U n ^Jn> J— ni Tn(J) ) Un y^nj J— m T n (J) 

Un ^(Snj a n ), J— ni T n (S, a)^ — U n ^(S n , a n ), J_ n , T n (S, [a n , a_ n ]) 
~\~ U n I (S n , a n ), J_ n , T n (S, [a n , a_ n ]) J — U n I (S n , a n ), J_ n , T n (S, a 



(42) 
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We show that the following identities are true 



Rmii^nj a n)? J— n) — -Rm((Sn> a n)) J— n)i V 1X1 ^ 71 

T,, n (S,a) = T,, n (S, [a n ,a_ n ]), V q G Q. 



(43) 
(44) 



This set of equations implies that if user n unilaterally switches from BS a n to a n but keeps its covariance matrix 
unchanged, then all other users' transmission rates as well as the prices charged for user n remain the same. 

Identity (l43l is straightforward as the interference at user receiver caused by user n is unchanged as long as 
user n's transmit covariance matrix remains to be S n . 

To verify (|44l . we first recall that the prices are defined as follows 



1'qr I n(S, a) 



E 



df m {Rn 



meJV g (a)\« 



OIL, 



R m =R m (S) 



H^„C m 1 (S_ m ) (l]i q + H g ^ m S m H^ m C m 1 (S_ r 



x H g .„ i S„ l H^ m C m 1 (S_ m )H Q 



(45) 



Take any BS q that satisfies q / a„ and q / a re . Clearly we have {m : m G jV 9 (a),m ^ n} = {m : m G 
A/"g([a n , a_ n ]), m ^ n}. For BS g = a„ or g = a n , although user n has changed its association, the other users' 
associations remain the same. That is, we again have {m : m G J\f q (eL),m ^ n} = {m : m G A/" g ([a„, a_ n ]), m / 
n}. Combining the above two observations with the fact that the transmit covariances S of all the users remain the 
same, we conclude that T g n (S,a) = T gjn (S, [a n ,a_ n ]), for all q G Q. This proves (1441 . 
Now using (l44l) . the first difference in (l42l i becomes 



U n ^(Snj a n)j J— ri) T n (S, a)J J7 n ^(S n , a n ), J_ n , T n (S, [a n , a_ n ]) 

C^n ( (Sn) a n)) J— id T n ,(S, [a n ,, a_ n ]) I — U n I (S n , a n ), J_ n , T n (S, [a n , a_ n ]) 



(i) 

— f ( (Sri j a n) ) J — n ) — / ( (Sri) a n) j J— n 



(46) 

where the inequality (i) is due to (|28l ), which states that for a fixed system association profile ([a n ,a_ n ] in this 
case), the user n's increase of utility induced by unilateral change of its transmit covariance is upper bounded by 
the increase of the system sum utility. The second difference in (l42l ) becomes 



^(S rl , a„), J— ri; T„(S, [a n , a_ n ])^ U n ^(S„, a„), J_ n , T„(S, a)^ 



= /„ (R n {{S n ,a n ),3- n ) -Tr T„(S, [a„, a_„])S„ - /„ i?„((S„,a„), J_„) +Tr T n (S,a)S„ 



— fn [ Rn ((Sri) a n ), J_n) ) /ri ( -Rri ((S r i, a n ), J , t 



(ii) 



iV 



^ fm ( -^nt((Srrt) a m ), J — m) ) /m ( ( (S m , a m ) , J_ r n ) J — / f (S n , a n ), J_ n j / ( (S n , a„), J_ 



(47) 
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where the step (i) follows from (l44l ) and step (ii) is due to d43l ). Combining (l42l ). d46l) and d471) , we have that 

U n f (Snj a n) i J— nj T n (S, a)^ — U n ^(S n , a n ), J— n, T n (S, a)^ 

< / ((S„,a n ), J_ n ) - / ((Sn.Sn), J- re ) • (48) 

This proves the claim. ■ 

Appendix E 
Proof of Theorem @] 

By the generalized potential property stated in Proposition^ the sequence {/(S* , a')} is monotonically increasing 
and converges. Let us denote the limit of this sequence as /*. 

Let A be the set of association profiles that appear infinitely often in the sequence {a*}. Take any a G A, define 
the subsequence {a* fe } of {a*} that satisfies 

a tk = a, and a m / a, V m G (t k ,t k+1 ). 

Clearly the sequence {f(S tk , & tk )}k is also increasing and converges to /*. Let S* be a limit point of {S* fc }. Take a 
further subsequence {a* fc } of {a tfc } such that lim^oo S tk = S*. Due to the fact that T n (-) is a continuous function 
in S, the subsequence {T^(S* fc )} must be convergent for all n. Let T* = T n (S*) for all n. We wish to show that 
(S*, a, T*) is a NE of the game Q J . The desired result will be shown in two steps. 

51) C7 n (S*,a,T£) > C/ n ([S n ,S* n ],[a n ,a_ n ],T*), V S n G T n Vn G M. 

52) £7 n (S*,a,T£) > [7 n ([S n ,S*„ n ], [a n ,a„ n ],X*), V a n ^ a n , V S n G F n . Vn G M. 

Without loss of generality (by possibly restricting to a further subsequence), we can assume that at time instance 
+ 1, it is user l's turn to act. 

Step 1) Let {S\) tk be the (unique) solution to the problem maxs ie j^ C7i(Si, S^, a tk , T* fe ). To show SI), it is 
sufficient to show that Hindoo ||(S*)^ — S* fc || = 0. Due to the strict concavity of E/i(Si, S_i, a, Ti) in Si, and 
use the definition of (S|)* fe , we have: 

Ui ((Si) t \S%SL t ",1^) > Ui (Si + t((S^* fe - Si^S^a'*,!?**) , V t G [0, l), Si + (S\) tk . (49) 

We will show that lim^^S*)** - S** = 0. The proof is along the lines of that of (39j Proposition 2.7.1], but 
with some important modifications, due to the lack of concavity/convexity of the function /(•) with Si. Suppose 
{(S^) tfc — S\ k } does not converge to 0. Let j tk = ||(S^)* & — S\ k ||, then by possibly restricting to a further subsequence 
of {t k }, we can find a 7 > such that j tk > 7, V k. Let V** = ((S^) t ' i - S\ k )/ r y tk , which is equivalent to 
(S^) ifc = S\ k + j tk V\ k . Clearly, ||V^|| = 1, and by possibly restricting to a further subsequence, we assume that 
V^ fc converges to Vi. 

Let us fixed some e G (0, 1). We must have < < r y tk . So S* fe + ejV^ lies on the line segment joining S* fc 
and S* fc +7 ifc V^ fc = (S*)**. Since (Si fc+1 ,a* k+1 ) is the solution to user l's utility optimization problem at time 
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t k + 1, we have (recall that J*^ = (S^a^)) 

^((Si fc+1 ,a^ +1 ),J^,T^) 

> Ui ( (SI ) tk , S\ , a tk , T\ k ) = Ui (S\ k + 7 * fc , , a tk , T\ k ) 

> U x (S ifc + e 7 V* fc , S\ , a tk , T[ k ) > Ui (S[ k , S\ , a tk , T[ k ) . 

where the last two steps follow from the concavity of U (cf. d49l)). Combining the above inequality with (|48T ). we 
have 

f((Sf+\a{* +1 ),J t * 1 ) - f((S{\4%J t * 1 ) 

> U 1 ((S t 1 h+1 ,a\ k+1 ), J%,T[*) - ^((S^), J*_? 1} T*») 

> ^i(S' fc +e7V' fc ,S^ 1 ,a <fc ,T f 1 fc ) - ^((S^.a^),^,^*) > 0. (50) 
Due to the fact that the sequence {/(S',a*)} converges to /*, we can take the limit of (l50l . and obtain 

0>^i(S; + €7V > Sl 1> a,Tl)-E7 1 (S*,a l T;)>0. (51) 

From the assumption, 7 > 0, and e G (0, 1), we have 67V 7^ 0. This contradicts the fact that for fixed a, 
J7i(Si, Sl 1; a, T*) has a unique maximizer (which can be seen by setting t = in (l49l)). We conclude that 
(S^)* fc — S** converges to 0. Due to the fact that lim^oo S tk = S*, we have (Si)' fc converges to S\. This implies 

E7i([S;,S!L 1 ] > a,Tl) > tfi([Si > Sl 1 ],a,T;) ) VS l£ Ji. (52) 

Step 2) From (l48i we have that 

/((S« k+ \a** +1 ), J*) - /((S**,a**), J*) 

> ^((S^^.a^ 1 ), J^.Ti*) - ^((S^a**), J*_\,Tf ) > 0. 

Utilizing the above relationship as well as the fact that the sequence {/(S*, a*)} converges to /*, we have the 
following limiting arguments (notice the fact that a' fc = a, V k) 

lim U 1 {S tk+1 ,a tk+l ,T tk )= lim £7i(S <fc , a* fc , T' fc ) = f7i(S*, a,T*). 

k— >oo k— >oc 

From the definition of (S^ fc+1 , a' fc+1 ), we must have 

U l {S tk+ \a tk+ \T tk ) > ^([Si.SL*!], [Si.a^J.Ti*), Va^ a\ k , V Si E Ji. 
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Taking limit of both sides, and notice the fact that a* fc = a, V k, we have 

lim U 1 (S tk+1 ,a tk+1 ,T t 1 k ) = E7i(S*,a,Tj) 

k— >oo 

> lim ^([Si.S^lJa^a^l.Ti*) 

k— >oo 

= C7i([S 1 ,SI 1 ],[a 1 ,a_ 1 ],T^), V ai ^ ai, VSiGJi- 

This says 

[7i(S*,a,Tt)>C7 1 ([S 1 ,S* 1 ],[a 1 ,a_ 1 ],TD, V a x / a ls VSjeJi. (53) 

Combining d52l ) and d53l , we have 

U l (S*,a,Tj) > U l ([Si.S*!], [ai,a_i],Tl) , V Si € Ji^ G Q. (54) 

Enumerating the above steps for all n G A/", we have that (1541 is true for every user, thus (S*,a, T*) is a NE of 
game Q J . 
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